DRAFT for lUTAM symposium 1999 



Version of February 8, 2008 
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The problem of anomalous scaling in passive scalar advection, especially with (5-correlated veloc- 
ity field (the Kraichnan model) has attracted a lot of interest since the exponents can be computed 
analytically in certain limiting cases. In this paper we focus, rather than on the evaluation of the 
exponents, on elucidating the physical mechanism responsible for the anomaly. We show that the 
- -- - anomalous exponents stem from the Lagrangian dynamics of shapes which characterize configu- 

, rations of n points in space. Using the shape-to-shape transition probability, we define an operator 

■ whose eigenvalues determine the anomalous exponents for all n, in all the sectors of the SO(3) 

' symmetry group. 
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en ■ I. INTRODUCTION 

I In the lecture at the lUTAM symposium the work of our group on the consequences of anisotropy on the universal 
. statistics of turbulence has been reviewed. This material is available in print, and the interested reader can find it in 
10^ ■ ^ short review is available in the proceedings of "Dynamics Days Asia" . In this paper we review some recent 
_ ^ , work aimed at understanding the physical mechanism responsible for the anomalous exponents that characterize the 
' statistics of passive scalars advected by turbulent velocity fields. We will consider isotropic advecting velocity fields, 
but will allow anisotropy in the forcing of the passive scalar. In such case the statistical objects like structure functions 
and correlation functions are not isotropic. Instead, they are composed of an isotropic and non-isotropic parts. We 
overcome this complication by characterizing these functions in terms of the S0(3) irreducible representations. Any 
^ r such function can be written as a linear combination of parts which belong to a given irreducible representation of 
SO (3). We will show that each part is characterized by a set of universal scaling exponents. The weight of each part 
however will turn out to be non-universal, set by the boundary conditions. 

The SO (3) classification will appear to be natural once we focus on the physics of Lagrangian trajectories in the 
flow. We will see that one can offer a satisfactory understanding of the physics of anomalous scaling by connecting the 
' the statistics of the passive scalar to the Lagrangian trajectories. This connection provides a very clear understanding 
[ of the physical origin of the anomalous exponents, relating them to the dynamics in the space of shapes of groups of 
Lagrangian particles. 
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II. THE KRAICHNAN MODEL OF PASSIVE SCALAR ADVECTION 

The model of passive scalar advection with rapidly decorrelating velocity field was introduced by R.H. Kraichnan 
already in 1968. In recent years |p|-p^ it was shown to be a fruitful case model for understanding multiscaling in 
the statistical description of turbulent fields. The basic dynamical equation in this model is for a scalar field T(r,t) 
advected by a random velocity field u{r,t): 

[dt-KoV^ + u{r,t)-V]T{r,t) = f{r,t) . (2.1) 

In this equation f{r,t) is the forcing and kq is the molecular diffusivity. In Kraichnan's model the advecting field 
u{r, t) as well as the forcing field /(r, t) are taken to be Gaussian, time and space homogeneous, and delta-correlated 
in time: 

((u"(r, t) - ^."(r', t)){uP{r, t') - uP{r', t'))) ^ = h^^{r - r')S{t - t') , (2.2) 
where the "eddy-diffusivity" tensor h°'^{r) is defined by 
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52„(i?)cxi?^-cx[^2(i?)]"(|) . (2.6) 



Here 77 and A are the inner and outer scale for the velocity fields, and the coefficients are chosen such that dah"^ = 0. 
The averaging (. . .)„ is done with respect to the realizations of the velocity field. 
The forcing / is also taken white in time and Gaussian: 

{f{r,t)f{r'X))^ = E{r-r')5{t-t') . (2.4) 

Here the average is done with respect to realizations of the forcing. The forcing is taken to act only on the large scales, 
of the order of L (with a compact support in Fourier space). This means that the function S(r) is nearly constant for 
r <ti L but is decaying rapidly for r > L. 

iFvom the point of view of the statistical theory one is interested mostly in the scaling exponents characterizing 
the structure functions 

^2„(ri,r2) = ((T(ri,t)-r(r2,t))2")^^^ . (2.5) 
For isotropic forcing one expects S2n to depend only on the distance R= \ri — r2| such that in the scaling regime 

In this equation we introduced the "normal" (^^2) and the anomalous {6271} parts of the scaling exponents = 
n(^2 — ^in- The first part can be obtained from dimensional considerations, but the anomalous part cannot be guessed 
from simple arguments. 

When the forcing is anisotropic, the structure functions depend on the vector distance R = r\ — r2- In this case 
we can represent them in terms of spherical harmonics, 

52„(ii)=E 

a£,m(-R)y«,m(-R) , (2.7) 

l.m 

where i? = R/ R. This is a case in which the statistical object is a scalar function of one vector, and the appropriate 
irreducible representation of the SO (3) symmetry group are obvious. We are going to explain in the next section that 

the coefiicients ai^rn{R) are expected to scale with a universal leading scaling exponent . The exponent will turn 
out to be £ dependent but not m dependent. 

Theoretically it is natural to consider correlation functions rather than structure functions. The 2n-order correlation 
functions are defined as 

F2„(ri, . . . , r2n) ^ (r(ri)r(r2) . . . T{r2n))^j ■ (2.8) 

For separations rij the correlation functions converges to (T^")^, whereas for r^- L decorrelation leads to 
convergence to (r)^ ^. For all rij k, 0{r) <C L one expects a behaviour according to 

F2„(ri,..., r2„) = L"(2-?)(co + ... + Cfe(r/i)?^"F2„(ri,...,f2„) + ---) , (2.9) 

where F2n is a scaling function depending on which denote a set of dimensionless coordinates describing the 
configuration of the 2n points. The exponents and scaling functions are expected to be universal, but not the c 
coefiicients, which depend on the details of forcing. 

It has been shown |^ | that the anomalous exponents C,2n can be obtained by solving for the zero modes of the exact 
differential equations which are satisfied by F2n- The equations for the zero modes read 

[-K^V,2+^2„].F2„(ri,r2,...,r2„) = . (2.10) 

i 

The operator B2n = Y^^>j ^ij' ^^'^ ^ij ^'^'^ defined by 

= B[r,,r,) = /i"^(r, - r,)9V9rf 9rf . (2.11) 
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III. LAGRANGIAN TRAJECTORIES, CORRELATION FUNCTIONS AND SHAPE DYNAMICS 



An elegant approach to the correlation functions is furnished by Lagrangian dynamics ||7|,|15H18[ . In this formalism 
one recognizes that the actual value of the scalar at position r at time t is determined by the action of the forcing 
along the Lagrangian trajectory from t = —oo to t: 



with the trajectory r{t) obeying 



r(ro,to)= / dt{f{r{t),t))^ , (3.1) 



rih) = ro , 

dtr{t) = u{r{t), t) + V2^r]{t) , (3.2) 



and ?7 is a vector of zero-mean independent Gaussian white random variables, (ri"{t)ri^ (t') ) = 5"^5{t — t'). With this 
in mind, we can rewrite -Pjn by substituting each factor oiT{ri) by its representation ( |3.l| ). Performing the averages 
over the random forces, we end up with 

F2n{ri,...,r2n,ta) = ( f dti ■ ■ ■ dt„[E{ri{ti) ~ r2{ti)) ■ ■ ■ (3.3) 



x^(r2„-i(t„) - r2„(t„)) + permutations ) , (3.4) 

To understand the averaging procedure recall that each of the trajectories obeys an equation of the form ([3.2|), 
where u as well as {ri^}^2i are independent stochastic variables whose correlations are given above. Alternatively, we 
refer the reader to section II of , where the above analysis is carried out in detail. 

In considering Lagrangian trajectories of groups of particles, we should note that every initial configuration is 
characterized by a center of mass, say R, a scale s (say the radius of gyration of the cluster of particles) and a shape 
Z. In "shape" we mean here all the degrees of freedom other than the scale and R: as many angles as are needed to 
fully determine a shape, in addition to the Euler angles that fix the shape orientation with respect to a chosen frame 
of coordinates. Thus a group of 2n positions {r^} will be sometime denoted below as {R, s, Z}. 

One component in the evolution of an initial configuration is a rescaling of all the distances which increase on the 
average like t^/''^; this rescaling is analogous to Richardson diffusion. The exponent C2 which determines the scale 
increase is also the characteristic exponent of the second order structure function S. This has been related to the 



exponent ^ of (2.3) according to C2 = 2 — ^. After factoring out this overall expansion we are left with a normalized 
'shape'. It is the evolution of this shape that determines the anomalous exponents. 

Consider a final shape Zq with an overall scale sq which is realized ai t — 0. This shape has evolved during negative 
times. We fix a scale s > sq and examine the shape when the configuration reaches the scale s for the last time before 
reaching the scale sq. Since the trajectories are random the shape Z which is realized at this time is taken from a 
distribution ^{Z; Zq, s — > sq). As long as the advecting velocity field is scale invariant, this distribution can depend 
only on the ratio s/sq. 

Next, we use the shape-to-shape transition probability to define an operator 7(s/so) on the space of functions '^(Z) 
according to 

ms/so)mZo) = J dZ-f[Z- Zo, s ^ so)*(Z) (3.5) 

We will be interested in the eigenfunction and eigenvalues of this operator. This operator has two important properties. 
First, for an isotropic statistics of the velocity field the operator is isotropic. This means that this operator commutes 
with all rotation operators on the space of functions "^{Z). In other words, if Oa is the rotation operator that takes 
the function 'i!{Z) to the new function \E'(A^^Z), then 

Oa7 = • (3.6) 

This property follows from the obvious symmetry of the Kernel ^{Z; Zq, s sq) to rotating Z and Zq simultaneously. 
Accordingly the eigenfunctions of 7 can be classified according to the irreducible representations of SO (3) symmetry 
group. We will denote these eigenfunctions as Bqim{Z). Here ^ = 0, 1, 2, . . ., m ~ —£+!,...£ and q stands for a 
running index if there is more than one representation with the same £,m. The fact that the BqgmiZ) are classified 
according to the irreducible representations of SO (3) in manifested in the action of the rotation operators upon them: 
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qim' 



(3.7) 



where DUr„ii^) is the SO (3) £ x £ irreducible matrix representation. 

The second important property of 7 follows from the (S-correlation in time of the velocity field. Physically this means 
that the future trajectories of n particles are statistically independent of their trajectories in the past. Mathematically, 
it implies for the kernel that 



7(Z; Zo, s ^ So) = J dZij{Z;Zi,s ~* si)j{Zi;Zo,si So) , s > si > sq 
and in turn, for the operator, that 

7(s/so) = 7(s/si)7(si/so) . 



(3.8) 



(3.9) 



Accordingly, by a successive application of 7(s/so) to an arbitrary eigenfunction, we get that the eigenvalues of 7 



have to be of the form aq,i = (s/so)''^" 



{±)C2''/' Bge,niZo) = I dZjiZ;Zo,s^ so)Sg,„(Z) 
So 



(3.10) 



Notice that the eigenvalues are not a function of m. This follows from Schur's lemmas ||T^, but can be also explained 
from the fact that the rotation operator mixes the different m's (3.7): Take an eigenfunction Bq£,n{Z), and act on it 
once with the operator Oaj^s/sq) and once with the operator ■j{s/so)0\. By virtue of (3^) we should get that same 
result, but this is only possible if all the eigenfunctions with the sa me £ and the same q share the same eigenvalue. 

To proceed we want to introduce into the averaging process in (3^) by averaging over Lagrangian trajectories of 
the 2n particles. This will allow us to connect the shape dynamics to the statistical objects. To this aim consider any 
set of Lagrangian trajectories that started a,t t — — cx) and end up at time t = in a configuration characterized by a 
scale So and center of mass Rq = 0. A full measure of these have evolved through the scale L or larger. Accordingly 
they must have passed, during their evolution from time t = —00 through a configuration of scale s > sq at least 
once. Denote now 



M2n(i, s So, Zo)dtdRdZ 



(3.11) 



as the probability that this set of 2n trajectories crossed the scale s for the last time before reaching sq, Zq, between 
t and t + dt, with a center of mass between R an d R + dR and with a shape between Z and Z + dZ. 
In terms of this probability we can rewrite Eq.(3.4) (displaying, for clarity, Rq = and t = 0) as 



F2„(-Ro = 0,so,Zo,t = 0) = J dZ J dt j dRn2n{t,R,Z;s^ sa,Za) 

E{ri{ti) - r2{ti)) ■ ■ ■ S(r2„-i(t„) - r2„(t„)) + perms (s; R, Z,t) 



dti ■ ■ ■ dt^ 



(3.12) 



The meaning of the conditional averaging is an averaging over all the realizations of the velocity field and the random 
rjj for which Lagrangian trajectories that ended up at time t = in ii = 0, so, Zq passed through R, s, Z at time t. 

Next, the time integrations in Eg. ( 3. 12 ) are split to the interval [— oo,i] and [i,0] giving rise to 2" different contri- 
butions: 



dti ■ 



dtr, 



dti 



dt:,. 



dtn + 



(3.13) 



Consider first the contribution with n integrals in the domain [— 00, t]. It follows from the delta-correlation in time of 
the velocity field, that we can write 

dti--- dtn S(ri(ti) - r2{ti)) - - - S(r2„_i(t„) - r2„(t„)) + perms (s; R, Z, t) 



dti--- dtn ^(riiti) - r2{ti)) - - - 5(r2„-i(t„) - r2„(t„)) -t- perms 



= F2a{R,S,Z,t) = F2nis,Z) . 



(3.14) 
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The last equality follows from translational invariance in space-time. Accordingly the contribution with n integrals 
in the domain [—00, t] can be written as 



j dZF2n{s,Z) J dtJdR fi2n(t,R,Z;s^ so,Zo) . (3.15) 
We identify the shape-to-shape transition probability: 

7(Z; Zo, s ^ So) = / dt dR fj.2n{t, R, Z; s ^ sq, Zq) . (3.16) 

J —oo J 

Finally, putting all this added wisdom back in Eg. ( 3.12 ) we end up with 

F2n{so, Zo)=I + j dZ-f{Z; Zo,s ^ so)F2„(s, Z) . (3.17) 

Here / represents all the contributions with one or more time integrals in the domain [t, 0] . The key point now is 
that only the term with n integrals in the domain [—00, t] contains information about the evolution of 2n Lagrangian 
trajectories that probed the forcing scale L. Accordingly, the term denoted by / cannot contain information about 
the leading anomalous scaling expo nent belonging to -Pjn, but only of lower order exponents. The anomalous scaling 
dependence of the LHS of Eq.( |3.17| ) has to cancel against the integral containing without the intervention of /. 
Representing now 

F2n{s,Z) = ^ ag^£„i(s)Bq£„(Z) , 

qlm 

I — '^^IqtmBqtmiZa) (3.18) 

qlm 



and substituting on both sides of Eq.( 3.17 ) and using Eq.( 3.10 ) we find, due to the linear independence of the 
eigenfunctions Bqhn 

aqXm{sa) = Iqtm + { — ] aqMnis) (3.19) 

To leading order the contribution of Iqe,m is neglected, leading to the conclusion that the spectrum of anomalous 
exponents of the correlation functions is determined by the eigenvalues of the shape-to- shape transition probability 
operator. Calculations show that the leading exponent in the isotropic sector is always smaller than the leading 
exponents in all other sectors. This gap between the leading exponent in the isotropic sector to the rest of the 
exponents determines the rate of decay of anisotropy upon decreasing the scale of observation. 



IV. CONCLUDING REMARKS 



The derivation presented above has used explicitly the properties of the advecting field, in particular the (5-correlation 
in time. Accordingly it cannot be immediately generalized to more generic situations in which there exist time 
correlations. Nevertheless we find it pleasing that at least in the present case we can trace the physical origin of 
the exponents anomaly, and connect it to the underlying dynamics. In more generic cases the mechanisms may be 
more complicated, but one should still keep the lesson in mind - higher order correlation functions depend on many 
coordinates, and these define a configuration in space. The scaling properties of such functions may very well depend 
on how such configurations are reached by the dynamics. Focusing on static objects like structure functions of one 
variable may be insufficient for the understanding of the physics of anomalous scaling. 
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